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X ■ Abstract The aim of this article is twofold. First we examine from a new angle the 

■ question of the recovery of time in quantum cosmology. We construct Green functions 

for matter fields from the solutions of the Wheeler De Witt equation. For simplicity we 
work in a mini-superspace context. By evaluating these Green functions in a first order 
development of the energy increment induced by matrix elements of field operators, 
we show that the background geometry is the solution of Einstein equations driven 
by the mean matter energy and that it is this background which determines the time 
lapses separating the field operators. Then, by studying higher order corrections, we 
clarify the nature of the small dimensionless parameters which guarantee the validity 
of the approximations used. In this respect, we show that the formal expansion in 
the inverse Planck mass which is sometime presented as the "standard procedure" is 
illegitimate. Secondly, by the present analysis of Green functions, we prepare the study 
of quantum matter transitions in quantum cosmology. In a next article, we show that 
the time parametrization of transition amplitudes appears for the same reasons that it 
appeared in this article. This proves that the background is dynamically determined by 
the transition under examination. 
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1 Introduction 



In General Relativity, the invariance under reparametrizations in time and in space 
leads to a set of local constraints. These impose that the total four-momentum of 
matter + gravity vanishes identically. Therefore, in contrast to the usual classical and 
quantum mechanics, there is no preferred reference frame to specify when and where 
events occur. This absence of external reference frame can also be conceived as arising 
from the dynamical character of space-time. Hence, since space-time is coupled to all 
matter transitions, it is coupled to the particular event we want to localize. Therefore, in 
a weak perturbation limit, the matter-gravity correlations resulting from the constraints 
should deliver the historical description of events performed in the background geometry. 
Then two questions arise: 

How to extract from the total wave function of the universe the properties of the 
background geometry in which time lapses are defined ? 

Under what circumstances is the extraction of a non-dynamical background a legiti- 
mate approximation ? 

Both questions have received attention and we refer to [IH @ H] for reviews in which 
various attempts are compared. We recall that it has been shown in how, 
starting from the solutions of the Wheeler-De Witt equation, the semi-classical character 
of gravity leads to the Schroedinger equation for matter fields. This treatment has 
been pursued in @] 110[ IHI and both higher order corrections and the criteria for 



having reliable approximations have been discussed. However there is no consensus 
on the manner to identify the background. In refs. 0|ll|[T^[|l^, the development 
in a power series of 1/G (Newton's constant) leads inevitably to backgrounds which 
are empty solutions of Einstein equations. Instead, in refs. proposed 
that the background solution must contain the mean matter energy as a source, in a 
manner similar to what happens in a Born-Oppenheimer treatment: in the adiabatic 
approximation, the protons feel the mean Coulomb potential delivered by the electrons. 
In this paper, we settle this debate and prove, in conformity with 0[§]|]1O[) that the 



background solution is the Einstein solution driven by the mean matter content. To 
reach this conclusion we proceed as follows. Instead of focusing on the wave function 
of the universe itself, in which matter and gravity are treated on a same footing, we 
center our attention on matter matrix elements evaluated at some given gravitational 
configuration. In particular, we analyze matter Green functions wherein the two inser- 
tions of the field operator are performed at two different radii of the universe. This 
un-symmetrical treatment in which a restricted number of matter quanta are involved 
determines unambiguously what is the background geometry and what are the develop- 
ments which lead to the localization of quantum matter effects. The small parameter 
controlling these developments is the relative change of matter energy induced by the 
insertion of the field operators appearing in the Green function itself. 

This is as it should be: given a certain matrix element, the validity of the background 
field approximation puts constraints on the sources of the background field, see Chapter 



14 in WM. Therefore the Planck mass has nothing to do with its validity nor does it 



govern its corrections. The corrections to the background field approximation are of 
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two kinds. One has the corrections to the WKB approximation for the propagation of 



gravity, see [Tl|, and the nonhnear dependence of the geometry in the energy change. 
At this point, it is relevant to notice the close analogy between the present development 
leading to the concept of time and the procedure by which a big isolated system furnishes 
a temperature to a little system contained in it. Both equilibrium concepts, i.e. time 
and inverse temperature, are obtained through a first order change in the energy and 
determined by the properties of the big system. Moreover, the second order corrections, 
which control the validity of these equilibrium notions, have the same origin: the big 
system is finite and has therefore a finite specific heat. 

Very important also is that only a small part of the corrections to the WKB approxi- 
mation affects the value of the Green functions. Indeed, only that part which is sensitive 
to the matter energy change modifies the Green functions. This has two important 
consequences. First, for matter dominated universes, the modifications of the Green 
functions due to these WKB corrections are irrelevant since they are one higher order 
than modifications induced by the non-linearities. Secondly, the validity of the WKB 
approximation is much less restrictive that the one applied on the full wave function of 
the universe. 

Finally, we clarify the role of the spread of the matter energy in the universe. In 



former works, see e.g. ||T0[1T^, the problem of the recovery of time was considered 



together with the validity of representing gravity by a single background. In this work, we 
emphasize that these two problems should be considered separately. Indeed we show that 
each component of the wave function defines its own time lapse. Therefore large spreads 
are perfectly acceptable. Simply the notion of a single background defining a single time 
valid for all configurations looses sense. But this has no physical consequences, since 
remote configurations do not interfere. 

In the next paper p^j, we define and evaluate the amplitudes of matter transitions in 
the framework of quantum cosmology. Upon evaluating these amplitudes, we shall see 
that a time parametrization is provided for the same reasons that time appeared in our 
Green functions, i.e. through a first order change of the energy change induced by the 
transition itself. We emphasize this point: Time appears because there is a transition, 
more precisely, time appears because there is an interaction hamiltonian which deter- 
mines the probability amplitude to have a transition. This legitimizes dynamically the 
analysis of the present paper which is purely kinematical. 

In this paper we proceed as follows. We consider two different matter fields, a massive 
one and a massless one conformally coupled. This allows us 

a. to reveal how, by using gauge independent quantities, proper time lapses char- 
acterize Green functions of the massive field whereas it is the conformal time which 
appears in the massless field's functions; 

b. to show how the appearance of these time lapses relies on a double development: 
a WKB approximation for the kernel of gravity together with a first order development 
in the change of energy induced by the field operator involved in the Green function; 

c. to show how the spectator field, not directly involved in the Green function, 
determines indirectly this function through its energy. Indeed, its energy contributes 
to the determination of the background geometry in terms of which time lapses are, in 
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turn, defined; and 

d. to analyze the corrections to these approximations and the physical conditions 
under which these corrections can be correctly discarded. 

All of these developments are done in the mini-superspace context, i.e. both the field 
configurations and the field operators are homogeneous and carry zero momentum. At 
the end of the article, we consider an extension of the formalism in which the momenta 
no longer vanish. 



2 Why study Green functions ? 

In this section we compare the time dependence of the Green functions of an harmonic 
oscillator to the time dependence of its wave function. The technical reason for compar- 
ing these well known features is to make the reader aware of the different roles that the 
total matter energy and the energy change induced by the insertions of the fields oper- 
ators will play upon considering, in Section 4, Green functions from the solutions of the 
Wheeler De Witt (WDW) equation. In that Section, we shall see that the total matter 
energy determines the background solution whereas the energy change characterizes the 
phase of the Green function. 

The physical justification for considering Green functions is the usual one: their study 
prepares the analysis of transitions amplitudes among matter constituents in quantum 
field theory. To study these transitions amplitudes in the framework of quantum cos- 
mology is precisely the subject of . 



In non-relativistic quantum mechanics, in the Schroedinger picture, the time depen- 
dence of the state of the system is governed by the following equation 

= e-^^(*-*o)/^|^o) (1) 

where |\l/o) is the wave function specified at t = to- 

When the system is in an eigenstate of the energy, this equation reduces to 

= e-^^(*-*°)/^|^o) (2) 

where E is the total energy of the system. For instance, for a harmonic oscillator in the 
n-th level, one has E = En = hu [n + 1/2). 

We now consider the time dependence of Green functions. This is usually done 
in the Heisenberg picture upon dealing with the time dependent Heisenberg operators 
g(t) = e^^ii-io)/hq^-iH(t-to)/n ^ p^^, instance, for an harmonic oscillator, whose state at 
t = to is l^o), the Green function is given by 

^,(^2,^1) = (^o|g(t2)g(ti)|^o) (3) 

When |\l/o) is the eigenstate |n) of energy En, one finds 

G^(t2, ti) = ^{n e*-(*2-*i) + {n + l) e-*"^*^-*^)) (4) 
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The differences between state propagation in eq. (0) and Green function time dependence 
are manifest upon computing eq. @ in tlie Scliroedinger picture: 

Gq{t2,ti) = iji\e^H{to-t2)/Ti q ^-iH(t2-ti)/h ^ ^-iH{ti-to)/n\^^\^ 

— (n\q g~*-f^(*2-ii)/?i q g-«-f^{ti-t2)/^|^^ 

= — 1?7,) e"*'-'^""^'^-"-*^"*^-'/^ + (r2| — |n) e~*(-^""'"^'^^*^*2-ti)/?i ^iE„(ti-t2)/h 

wliere we liave introduced tlie (Scliroedinger) time independent expectation value {n\aa^\n) 
n + 1. The total energy En drops out the phase of the Green function. Only ±a;, the 
energy changes induced by the q operators, survive. This will be extremely important 
in the quantum cosmology. In that context, the Green functions are delivered in the 
Schroedinger form, as in the third line of eq. (^). Furthermore, the two independent mat- 
ter energy scales hu, the (microscopic) change in energy induced by the operator g, and 
En = nhuj, the (macroscopic) energy in the universe, play completely different roles0. 
Indeed, it is by a development in 1/n that the background solution characterized by the 
total energy is isolated. Furthermore, higher order terms in 1/n determine the correc- 
tions induced by the fact that geometry is dynamical and can only be approximated by 
a single passive background. 

The other difference is that the Green function is independent of to, i.e. of the initial 
phase of the wave function and depends only in the lapse between t2 and ti. This 
local character will also play an important role in the WDW context since only local 
corrections will enter into the nonlinear response of gravity. 

However it should be pointed out that upon dealing with a superposition, |\E'o) = 
^nCn\n), a dependence in to and t2 + ti does appear. Indeed one obtains 



G';"(t2,ti) = S„|c„|2i- (n e*"(*^-*^) + (n + 1) e'^'^^*^-*^)) 

1-Cnc*n_2^n{n - 1) e— 



2uj 
+S„2 Re 



(6) 



The dependence on and t2 + ti emerges from the interfering terms relating |n) to 
|n ± 2). Therefore, when the system is described by a density matix which is diagonal 
in energy, i.e. p = E„p„|n)(n|, the interfering terms do not exist and one has 

G^(t2, h) = Tr[ p g(t2)g(ti)] = S„p„^ (n e*-(*^-*i) + (n + 1) e-'^^^'^-'^^) (7) 

Because of this simplicity, in Section 4, we shall first consider Green function arising 
from diagonal density matrices and only then discuss the new aspects brought in by the 
interfering terms. 

3 The cosmological model 

We recall the main features of the dynamics in mini-superspace, see e.g. We put 

special emphasis on the relationship between the matter energy content of a universe and 
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lapses of proper time evaluated in that universe. The reason for this emphasis is that 
the time parametrization reappears in quantum cosmology through this relationship. 

To make this as clear as possible, we choose the matter content in a such a way that 
its state is characterized by constants of motion. This allows to write the Hamilton- 
Jacobi action as a sum of actions and the wave function as a sum of products. Because 
of the simplicity of this model, the modifications introduced by the quantum character 
of gravity will be easily identified and evaluated. More elaborate situations will be 
consider afterwards. In Appendix A, we leave the mini-superspace restriction and in 
we introduce interactions which cause the "constants of motion" to evolve. 

We consider a cosmological model with both massive and massless matter fields. For 
the massive field, we impose that the mass M be large enough that its Compton wave 
length satisfies 

h a 

M « -a W 

where a = da/dt with t being the proper time of the universe, i.e. when the lapse satisfies 
= 1. We recall that, in mini-superspace, the metric element reads 

ds^ = -N^iOd^ + a^Od^n^ (9) 

where N{^) is the lapse and where d^Q^ is the constant line element of the homogeneous 
three surfaces. 

When eq. is satisfied, there is no pair production of massive quanta and the 
action of an homogeneous massive scalar field ip{t) is correctly approximated by (from 
now on we put h = 1) 

SM = \j dt _ mV') (10) 

where the field il){t) is related to the original scalar il){t) by il){t) = a^^'^{t)ilj(t). Indeed 
one verifies that the differences between the original field equation and the harmonic 
(adiabatic) equation for ilj{t) scale like (a/aM)^ << 1. Non adiabatic corrections, which 
govern pair creation amplitudes, will be considered in |2l|. 

For the massless conformally coupled field, see |T6|, its action is 

S. = lldv {{d,<Pr - k'<P') = 11^ {a^' - k^') (11) 

where drj = dt/a{t) defines the conformal time. It that gauge, N{ri) = a{r]). The field (j) 
is related to the original massless field (p hj (p = acf). 
The total energy of matter (when A^ = 1) is given by 

H^{a)=M\AM\'' + -\A^\^ (12) 

(X 

where Am,A^ are the classical amplitudes of the ip field and the (p field whose norms are 
conserved, i.e. are a independent. 

The action for gravity is given by the Einstein lagrangian. When restricted to the 
conformal mode in the proper time gauge is reduces to 



S, 



G 



1 f 

—— / dt (—a d^ — Ka — Aa^j (13) 
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where k is equal to 0, ± for flat, open or closed three surfaces, where A is a cosmological 
constant and where G is Newton's constant. Notice the minus sign of the kinetic term 
of a. This is a common feature of reparametrization invariant theories. In the present 
case, the total hamiltonian of gravity plus matter is constrained to vanish in order to 
guarantee the invariance under time reparametrizations. In terms of the momentum of 
gravity, tTq, this constraint is given by 

where Ha is related to a and to the proper time t by 

71"a = (15) 

The use of the momentum tt^ allows to write the constraint without any reference to 
the lapse A^, i.e. in a gauge independent manner. Thus, for a given matter content, in 



our case with Am, specified, tt is determined (up to a sign) by the constraint eq. ([1^) 



7r(a, Em, E^) = -G'^^na^ + Aa^ + 2Ga (^Em + ^ 



(16) 



where Em = MIAmI"^ is the energy of the massive field and E^ = A;|A^p the "conformal" 
energy of the massless field. The sign ambiguity is fixed by restricting ourselves to 
expanding universes, i.e. d > 0. Therefore, instead of determining tt in terms of a given 
a{t), eq. (^) should be considered as determining the proper time t in terms of n and 
a having specified the matter content. Indeed one has 

r"2 , a 



ra.2 

At{a2,ai;EM,E^) = - da- 

Jai (jr TTyCl, iLiM, ^'y) 

= -9em / da 7r(a, Em, E^) (17) 



Similarly, the lapse of conformal time Ar] is determined by 

r«2 1 



ra.2 

Ar]{a2,ai;EM,E^) = - da- 

Jai G 7r(a, Em, -c/7) 

/•a2 

= -dE.y / da 7r{a, Em, E^) (18) 

J ai 



In the second lines of eqs. (|T7|, |1^), we have expressed the fact that the proper time 
t is, by definition, the conjugate variable of the energy Em and the conformal time r] is 
the conjugate variable of the conformal energy i?^. In both cases, it is the gravitational 
part of the action, i.e. Sq = J da tTq, see below eq. (0), which is used to define the 
conjugate variable. By using the angle-action terminology, we can say that the proper 



time is the "angle" conjugated to the "action" Em, see e.g. [17| for a brief and clear 



presentation of these concepts. When using the solutions of the WDW equation to 
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evaluate Green functions in quantum cosmology, we shall see that it is through these 
equations that time lapses will be defined and determined, answering therefore point 
a. of the Introduction. Moreover, these Hamilton- Jacobi equations make manifest the 
analogy with the inverse temperature which is defined in a microcanonical ensemble as 
the derivative of the entropy with respect to the energy, see eq. (|5^. 

Since Am and are conserved, the Hamilton- Jacobi action is a sum of three actions 

S{a, ^, 0) = Scia; Em, E,) + ^^/-(V^; Am) + S';-^-(0; A,) (19) 

where the matter actions are the actions for harmonic oscillators of amplitude A and 
respective frequency M and k. The gravitational part of the action satisfies 

- G^{daSG{a)f + Kc? + ka^ + 2Ga [em + ^) = (20) 

Since Sq depends on a single variable only, it is entirely determined by the classical 
momentum, eq. (0), by 

ra.2 

SG{a;EM,E.y) = / da 7r{a,EM,E.y) (21) 

J ai 

Upon quantizing the system, the conservation of the energies Em and E.y implies 
that the wave function factorizes into a sum of products of waves 

^(a, V',0) = J:nj,j,n-,CnM,n.y'^{a;nM,n^){^\nM){(p\ny) (22) 

The matter states Ium) and |n^) are states of harmonic oscillators characterized by the 
(conserved) number of excitations um and n^. They keep their usual normalization and 
interpretation as providing the quantum state of the oscillators, i.e. the interpretation 
that one uses when working with quantum matter fields in a classical geometry ||16||. This 
remains true when the states are adiabatically distorted by gravity, see Appendix A. 

Then, in order for the coefficients c„j,,^^„^ to keep also their interpretation as the 
probability amplitude to find the matter state in the nM,n.y sector, a normalization of 
the gravitational part of the wave function, '^{a;nM,n.y) must be chosen. This wave 
satisfies the WDW equation at fixed um and 

G'^d'l + Ka^ + Aa^ + 2GaEm{nM, n.y, a) \I'(a; um, n.y) = (23) 

where Em{nM,n^,o) = M{nM + 1/2) + A;(?2^ + 1/2) /a is the eigenvalue of the matter 
hamiltonian, eq. ([T2|) . Our aim is not to provide a new interpretation to ^'(a; n^^, n^) 
nor to discuss the "normal ordering" problem, i.e. to define the operator d^. Instead 
our aim is to determine (i) how to normalize ^ in order to be able to interpret the 
CnM,n-^ amplitudes of probability (ii) under which conditions its properties lead to a 
localization of quantum matter transitions similar to the one used when one works in a 
background geometry and (iii) what are the origin and the meaning of the corrections 
to this approximate description. We refer to for discussions concerning possible 
interpretations of the solutions of the WDW equation. 
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We shall prove that the necessary condition to interpret the Cnj^^^n-y as amplitudes of 
probability is that the waves '^{a^nu^n^) have positive unit (conserved) wronskian, i.e. 

'^*{a;nM,nj) ida'^{a\nM-,n^) = 1 (24) 

for all nM,n^- When considering the solutions of the Klein-Gordon equation, this con- 
dition with t replacing a means, in the absence of pair production, i.e. when the WKB 
approximation is valid, that one deals with a single particle. Similarly, in the present 
case, it means that one works with a single universe. 

In order to prove that unit wronskian allows that interpretation, we shall not work 
with the functions \E'(a; um, n^) directly. Instead we shall work with a kernel of the form 
K = The reason for this choice is that matter matrix elements are governed by K 

and not \E'(a; riM; ^^7)- Indeed, both quantum matter field theory and the gravitational 
corrections induced by the dynamical character of a will be delivered by the properties of 
K and not those of \E'. Furthermore, the choice of unit wronskian is elegantly introduced 
and justified by gauge invariance when working with K. 

Indeed, in strict analogy with the Feynman boundary condition for relativistic parti- 
cle in the Schwinger formalism|T^, one introduces a "fifth" time, here it corresponds to 
a constant lapse N, and one writes the Schroedinger equation (first order) for the kernel 
K{a2, fli; um, ny, N) to propagate from 02 to ai in a time 

idN — G'^dl + Ko^ + Aa"^ + 2GaEm{nM,n^, a) K{a2, ai;nM,nY, N) = S{N)5{a2 — ai) 

(25) 

This dependent kernel satisfies the usual boundary condition 

^im^ K{a2, ai, um, n^', N) = 6{a2 — ai) (26) 

Then one chooses the integration contour of N. In strict analogy to the Feynman 
boundary condition, we integrate over positive only and define 

K{a2,ai;nM,n^) = dN K{a2, ai^nM^n^-, N) (27) 



where the factor of has been introduced for further convenience. 

We point out two properties of this kernel which arise from the second order character 
of the WDW equation and the integration over positive only. First, its normalization 
automatically leads to positivity unit wronskian, in the following sense 

lim {ida2 - ida^) K{a2,ai;nM,n^) = (28) 

a2— ai-+0+ 

Secondly, the "convolution" of two kernels satisfies 

K{a3,ai;nM,n^) = K{a3,a2;nM,n^) ida2K{a2,ai;nM,n^) (29) 

when 03 > 02 > ai . This relation coincides with the convolution in time of the Feynman 
Green function of a relativistic particle, and not with the convolution of non- relativistic 



9 



kernels. The reason for the discrepancy is the additional integration over the "fifth" time 
N which projects onto zero energy solutions. We shall see in the next Section how these 
two properties guarantee that quantum matter field theory in a geometrical background 
will re- appear in a well defined scheme of approximations. In this respect, there is no 
necessity to provide an interpretation to K{a2, ai). To determine how to normalize and 
to use it is sufficient. 

Finally we rewrite this kernel introducing the density matrix of the matter. The 
reason is to unify the notations for matter Green function together with propagation of 
gravity. In the case of eq. (^), the density matrices are trivial: pm = |?T-A/)(^Af| and 

= \n^){n^\, but the following expression is valid for any density matrix. We express 

K{a2,ai;pM,P-y) as 

K{a2,ai;pM,P'y) = Tr[ pMP^\ai){a2\ ] (30) 

where the symbol Tr means the trace over all matter + gravity configurations such that 
the WDW constraint is satisfied by integration over positive A^. 

To get some fiavor of what is the physical meaning of the "Feynman" boundary 
condition in the present case, it is appropriate to consider the the WKB approximation. 
In this approximation, the kernel is equal to 

exp (i J^^ 7i{a, Um, n^)da\ 
K{a2,ai;nM,n^) = — , =^= = '^WKBia2;nM,ny)'^^j^^{ai;nM,ny) 

2^7r(a2, Um, n^)vr(ai, Um, n^) 

(31) 

where the momentum 7r(a, , ra^) is given in eq. (|T6D , and where '^wKB{0']nM,n.y) is 
the WKB solution of eq. ( |25D with unit wronskian. Therefore in that approximation, one 
has monotonic expansion when 02 > ai if the classical momentum tt^ does not vanish. 
Then, "negative energy" solutions, with vr > 0, are not generated since there is no back- 



scattering in the WKB approximation. Upon considering matter interactions in |]20 
negative energy solutions might be produced. We refer to that paper for a discussion of 
this point. 

In Section 5, we shall evaluate the corrections to the WKB approximation. For the 
moment, we just recall that its validity requires, see e.g. [ff6|, that 



daiT{a,nM,n^) « tt {a,nM,n^) (32) 

The reason for which we recall this criterion is the following: we want to compare the 
implications of eq. (^) and eq. (|]). To this end we introduce the dimensionless 
quantity Ng, given by Ns{a) = 2GEm{nM,n^,ci)/ci which may be understood as the 
"Schwarzschild radius" of the matter content of the universe measured in the unit of a. 
For a closed universe, one has A^^ = 1 at its turning point. More generally, in (non-empty) 
cosmology, one is confronted with "huge" numbers, i.e. numbers whose logarithm are 
much bigger than one, as well as large numbers, i.e. numbers whose logarithm are bigger 
but close to one. What we want to stress is that Ns belongs to the second class. (We 
shall consider separately the de Sitter case.) Thus upon dealing with "huge" numbers, 
like Um, one can put, in a first approximation, Ns = 1 and trade therefore a for the 
total matter content of the universe. 
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To see how this works, consider a matter dominated universe with Ns{a) = 0{1) and 
with Ns{a) ^ Aa^. Then, the condition that the Hubble radius be much bigger than 
the Compton wave length, eq. (||), leads to a constraint on um, the number of massive 
quanta, 

"« gSJ = ' (^^' 

Indeed one has {a/af ~ GMum/o? ~ [GMum)''^ when Ns = 0{1). The last inequality 
in eq. (^) simply recalls that M must be much smaller than the Planck mass Hpianck 
in order to legitimize the neglection of the gravitational dressing of individual quanta. 
This requirement has nothing to do with the semi-classical of the propagation of the 
cosmological radius a. In this respect, it should be pointed out that the condition (p3D, 
is not compatible with the limit G — which is used in refs. 0|l2|||l^ to define the 



background solution. This is related to the fact that we shall find, in contradiction to 
these works, that the background solution is driven by the mean matter energy. 
Using the same substitutions, eq. (|32D, gives 



» —Jf- (34) 

When eq. (|33|) is satisfied, eq. (|32D is fulfilled. This simply means that when the Hubble 
radius a/a is much bigger than the Compton wave length, matter dominated universes 
behave semi-classically. We emphasize this point: the validity of the WKB approxima- 
tion for a puts constraints on the sources of gravity, i.e. on matter. Furthermore, this 
hierarchy of constraints on um determines the relative importance of the corrections on 
the various approximations that we shall perform in the next section. This will be made 
explicit in Section 5. 



4 Green functions from WDW solutions 

To construct Green functions from the solutions of the WDW equation, we shall proceed 
in two steps. First, we work in the background field approximation, i.e. we evaluate the 
Green functions of ip and (j) from their time dependence in a classical universe whose 
radius follows a = a{t). Secondly, we define an extension of these Green functions in 
the WDW framework and show under which conditions the quantum matter-gravity 
evolution leads to the expressions previously obtained. 

In a classical universe whose radius is given, once for all, by a = a{t), the Green 
functions of the field are defined in total analogy with eq. (^, by 

G4t2,h) = (xo|^(t2)^(ti)|Xo) (35) 

where |xo) is the state of the ip field. When the state of the field is described by a 
diagonal density matrix, one obtains 

= S„,,p.,,-i- {um e^*^(*-*^) + {um + 1) e-^*^(*-*^)) (36) 
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where the symbol Tr^ has been introduced to be used in concordance with eq. (pOl). It 
means "trace over the configurations" . The last equality follows from the fact that the 
if) field is an harmonic oscillator of frequency M for all cosmological histories described 
by a(t), see eqs. (|, |^). For the field instead, one gets 

Gl'it^.tr- k) = S„^p„,^ {n, e^'^fe-''^) + {n, + 1) e'^'^^''^-''^)) (37) 

where the lapse of conformal time is defined by r]2 — r]i = J^^ dt/a{t). 

We now define a matrix element which generalizes eq. (|3^) in the WDW framework. 
In this case as well, we have to specify the matter content of the universe, i.e. the weights 
PnM^Pn-,- However we shall see that the matter content will now play a double role, the 
former one appearing in eq. (|36|) and the one of eq. {\17\} which determines the time 
lapse function of 02 and ai. 



Using the notation presented in eqs. (|30| , |36D , we define ^^(ci2, ^i) to be 



^/(a2,ai) = Tr 



pMP-y \a'2)ida2i^{a2\ \ai)ida^i/j{ai 



(38) 



The operator is the Schroedinger operator = [dM + d\.j)/V2M. Then the operator 

\ci'2)'ida2'4'{0'2\ should be read at a = 02". It replaces exactly the Heisenberg operator 

ilj{t2) in eq. ( ^61) which is at t = ^2"- The occurrence of the operator 27ra2 — ^^02 
follows from the second order character of the propagation of a, see the discussion after 
eq. (|2^) . We shall see that this operator guarantees that the weight attributed to "to 



be at a = 02" does not depend on the momentum at 02, hence "to be at a = 02" 
is independent of the matter content of the universe^ i.e. independent of nM,n^- This 
independence is required in order to interpret the matrices p of eq. ( |55| ) as providing 
the probability to find the system characterized by nM,n^, being at a = 02. 
In terms of the kernel at fixed matter content, see eq. (0), one gets 



G^{a2,ai) — '^^UM^n^PnuPn-f 



^^A(a2,ai;nM - l,n^) H ^j^K[a2,ai;nM + 1,^7) 

ida2ida^K{a2,ai;nM,n^) (39) 



where the two terms come, as in eq. (||), from the products d^j^dM and the dMd^^ which 
decrease and increase respectively the number of massive quanta by one unit. 

The origin of this operator 27ra can also be understood as foUows. Green functions are densities 
with respect to t (or a) since both are continuous variables. This means that Green functions enter 
into physical rates of transitions integrated over t (or a). Therefore instead of considering G^{t2,ti) 
one might consider 0^^(12, ti)dt2dti directly. Then, the corresponding quantity in quantum cosmology 
becomes G'^^ {a2,ai)da2dai/TTa2'^ai and the extra factors of 7r~^ cancel with the previous ones. Upon 
considering interactions, it is the hamiltonian framework which specifies the weight accompanying Green 
functions in transition probabilities, and indeed one finds that unit wronskian is required ||20| . However 
extra factor of l/-\/2|7r^ will still be found at the end points of the amplitudes. Their amputation is 
mandatory in order to obtain amplitudes of probability, as in the reduction formula. 
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Since um and are strictly conserved, decomposes as 

where ^^^^'"^(a2, Oi) is the Green function in a universe with exactly um and quanta. 
This is our first result. 

We now proceed to a double development. First, we use the WKB expression for 
K{a2,ai;nM,n^), eq. (|3T|), and we obtain (compare with eq. (|])) 

g^ n«2,ai) = ^ ^- ^ ' 

+ ^^^+e ^''i (41) 

where n{a,nM,n^) is the classical momentum, solution of eq. (p!6D, and where the 
normalization factors A± are given by 

1 [7i{ai,nM,n^) + vr(ai,nM ± l,n^)] [vr(a2, nM,n^) + n{a2,nM ± 1,?^7)] . 

=•= ~ 7 1/2" j 

4 [7r(ai, njv/, n^) vr(a2, nM, n^) 7r(ai, nM ± 1, n^) 7r(a2, ± 1, ^^7)] 

Secondly, we develop the phase and the norm of to first order in the change of 
the matter energy (= M), that is to the zeroth order in I/um- Using Hamilton- Jacobi 
relation, the phase can be written as 

[7r(a, njv/, "-7) — 7r(a, riM — 1, %)]c?a = j rr{a,nM,ny)da + O^I/um) 

= -M At(a2,ai;nM,r^7) +C(l/nM)(43) 



by definition of At{a2,ai;nM,n^), c.f. eq. {\l7[j. To this order in I/hm, one finds that 
the norm A± = 1, for all nM,n^. This is the necessary condition mentioned above. The 
gravitational propagation must be normalized in such a way that it gives a unit weight 
to the matter matrix elements, at least to this order in I/hm- 
Collecting the results, one gets 

= G;*^(At(a2,ai; nM,n^)) (44) 

see eq. (pG]). 

Therefore, we have proven that to first order in the change in the matter energy 
and to the WKB approximation of the kernels for gravity, ^^*^'"^(a2, ai), defined in eq. 



(|38|) , is equal to the Green function G^*^(At(a2, ai; nM,n^)) evaluated in the universe 
characterized by um and n^. This is our second result. From eq. (|4T|), we see that the 
sole role of the product of the two normalized kernels K{a2,ai] nM,ny) is to provide 
the phase of the Green function C^" {At{a2,ai; nM,n^))- 

The double dependence in um, that we mentioned after eq. (|37|), is now manifest. 



As in eqs. (^, ^), there is the usual "harmonic" dependence of the norm in the number 
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of quanta under examination, i.e. those of the field, but there is also the parametric 
dependence of the lapse of proper time At(a2, cii; nM, n^) in the total number of quanta 
of all species that are present in the universe. For exactly the same reason, i.e. the dy- 
namical character of gravity, a double dependence of the temperature in the occupation 
numbers is also found upon considering statistical mechanics in the presence of gravity, 
see the discussion after eq. (25) in . 



We emphasize that the quantity appearing in the Green functions, i.e. MAt = 
dn J daiTa is gauge independent. Indeed, it gives the number of nodes of the Comp- 
ton frequency from 02 to ai in a universe filled with um and quanta. The 

reason which guarantees this gauge invariance is that (i) the "entries" of the kernels 
K{a2,ai; nM,n^) are the gauge independent quantities a, tim and n^, (ii) the Green 
functions are determined in terms of variations of K with respect to these quantities. 

We also emphasize that the decomposition in eq. (|40|) requires only the existence of 
conserved quantities, in our case tim and n^0. No other approximation is needed. Eq. 
(|40|) is essential in guaranteeing that the background geometry, from which time lapses 
are defined, is the one characterized by these quantum numbers. In this we confirm and 
generalize [|10| and disagree with [|1^ |]13| . Indeed there is no reason to develop 



eq. (^Sp around um = = 0. It is this procedure which is (implicitly) adopted in 
those works by taking the limit G ^ in eq. (p^). 

In this respect, it is very instructive to compare the parallel development which 
leads to a canonical distribution for a little system contained in a much bigger isolated 
ensemble. First, one realizes that the action of the cosmological expansion, eq. (PH) , 
acts like the entropy of a reservoir in delivering the time lapse. Indeed, a heat reservoir 
determines the (inverse) temperature through a first order energy change of its entropy, 
compare eq. (^) and eq. ( |5^ in Appendix B. Secondly, upon considering second order 
energy changes, the analogy is maintained and reinforced since the corrections to the 
background field approximation are similar to to the finite size effects of a big isolated 
system, see Section 5. Thirdly, the additive character of the energy implies that two 
different massive fields evolve with the same proper time. Indeed, to first order in the 
energy changes, the derivative of eq. ( |i3| ) would determine the same time lapse for both 
fields. This is exactly like the "zero-th law" of thermodynamics (i.e. the equality of 
the temperatures at equilibrium) which follows in a microcanonical ensemble from the 
additivity of the energy and the equipartition ansatz. 

Before proceeding to the evaluation of the corrections and exploring the correlated nature 
of both approximations used, three important remarks should be added. 

1. Had we studied the corresponding Green function for the massless conformally 
coupled field 0, we would have obtained 

^?J*''''"(a2, = \ (n^ e'^^nia^^a,, „M,n^) + + 1) g-ifcAr,(a2,ai;nM,n7)A (45) 



instead of eq. (H). This is because the gravitational phase in the WKB approximation 



^ Wc shall relax the constraint on the Hubble radius, eq. (||), in and consider pair creation of mas- 
sive quanta. Then exactly conserved quantities will be replaced by adiabatically conserved quantities. 
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is governed by the change of / nda with respect of n^: 

ra.2 

dn^ 7r{a,nM,n^)da = -k Ar]{a2,ai;nM,n^) (46) 

by definition of Ar]{a2,ai;nM,n^), see eq. (0). Thus, the response of gravity, that is 
the first order change of / da 7i{a, hm, n^) induced by the operators -0 or 0, furnishes (by 
definition) the correct time dependence through Hamilton- Jacobi equations. 

2. Eq. (]40| ) tells us that the Green function is given by the weighted sum of Green 
functions at fixed nM,n^. Therefore one has 

^^(02,01) = ^nM,n^ PnuPn-, G'^'''{At{a2,ai;nM,n./)) 

= p„„ G;*^(At(a2, ai; um, n^)) = G^"(At(a2, ai)) (47) 

where nM,n^ are the mean number of massive and massless quanta and where the 
mean lapse At{a2,ai) is equal to At{a2,ai; fiM^n^y)- In the second equality, we have 
replaced the parametric dependence of Q through At{a2, ai.nM.n^) in um and by 
a single "background" contribution controlled by the mean occupation numbers. This 
is physically relevant for well peaked distributions, i.e. for {{AnY) / (n^) <<c 1. In that 
case, the Green functions Q are correctly given in terms of the properties of the mean 
universe only. This mean universe is the solution of Einstein equations driven by the 
mean energy density [f7||P| [p!0[] (the so-called semi-classical equations). In our case, it is 
determined by, see eq. (0), 

= 

■— + — + A + - 
a^ a^ a 




2 + — + A + — Mum + 7^7 ) = (48) 



where a means dia, i.e. derivative with respect to the mean proper time. 

Instead, for widely spread distributions, it is no longer meaningful to define a mean 
universe governing mean lapses because remote matter field configurations will classically 
evolve differently and quantum mechanically will never interfere. However this is free 
of consequences since each contribution at fixed nM,n^ delivers a well defined time 
dependence to Green functions. We shall return to these "decoherence" aspects upon 
considering interactions, in [pil[]. 

3. We stress that eq. (|41| ) delivers the Green function in the Schroedinger picture, 
see eq. (^. We now understand that the double dependence on um and the parametric 
dependence in lead to non-linearities which prevent to define an Heisenberg operator 
ip{a) from H]\j only. On the contrary, the superposition principle still survives for the 
Schroedinger evolution of the full wave function, here given by K{a2,ai;nM,ny). 



5 The gravitational corrections to Green functions 

There are three types of corrections associated with the three approximations used. 
First there are the nonlinear terms in the change of the matter energy induced by the 
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operators tp, see the passage from eq. (p| ) to eq. (0). Then one has the corrections 
to the WKB approximation used when replacing the kernels in eq. (|39|) by their WKB 
expression given in eq. (0). Finally one has the corrections associated with the spread 
in Um and upon computing the Green function in the mean universe, see eq. (^7]) . 
The first two are of dynamical character while the third one depends of the choice of 
the state of matter, i.e. the spread in pm and p^. Our aim is to determine what are the 
dominant and hence relevant corrections. 

These three corrections make explicit the non-linearities induced by the dynamical 
character of gravity. In particular, the approximate nature of the background is displayed 
and its sensitivity on the quantum transition under examination is pointed out. 

We start by evaluating the non linear corrections in the energy change which arise 



in the phase and the norm of eq. (|T). To second order in the energy change, the phase 
of the first term gives 
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ai 
0.2 



7r(a, Um, n^) — iT{a,nM — 1, '"■7)] da 
aM I aM 



ai |_G7r(a, nM, ''^7) 
MAt(a2, ai; tt-m, '^-7) 



aiM 
G'n'^{ai,nM-, ny) 



da 



1 + At(- 



27r. 

71 



(49) 



Notice that we have developed the expression around um and not around the "mean" 
value riM — 1/2. We have developed around hm since we want to use the same background 
for both terms of eq. (|4l|) . 

In the third line we have replaced the integral of the correction by the value of the 
integrand at ai times the interval plus the linear dependence in At. This is valid if a/vr^ 
is slowly varying. To estimate these terms, one should specify the dominant matter 
content of the universe. When the universe is matter dominated, the first correction 
term is 

^ M < J_ (50) 

G'7r^(ai, n^) 2MnM + '^kn-^/ a 2nM 

It is independent of G, h and a. It depends only on the ratio of the matter change 
(= M) and the total energy in the universe. This is what guarantees the validity of the 
background field approximation: that the sources of gravity in the universe be heavy[| 
when compared to the transition upon examination. This is our third result. 

The origin and the meaning of this term are clear. It arises from our choice to 
develop eq. (^91) around um and not around the "mean" value um — 1/2. Indeed, 
had one developed around the "mean" value, it would not be present. Its meaning is 
therefore that each term of eq. (0) defines its own background. This dependence of the 



background in the quantum event under examination will be further discussed in [20 



As already pointed out in Section 4, there is a close correspondence between the 
present development of the gravitational action and the development of the entropy of 

^As pointed out by R. Brout upon discussing these results, the emergence of an inertial time dehvered 
by the matter content of the universe can be viewed as rephrasing (and legitimizing) Mach's principle 
in this quantum context. 
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a big isolated system around equilibrium. Compare eq. (|i9| ) and eq. (0) in Appendix 
B. We reemphasize that both first order and second order energy changes are in strict 
correspondence. Indeed, first order changes determine time and inverse temperature 
and second order changes determine the corrections to these background-equilibrium 
concepts. 

We emphasize that the Planck length does not appear in this correction and we point 
out that this correction has the same origin as the recoil effect of a relativistic particle 
affected by a transition |]19[. In that case the correction terniQ is uj/2M where uj is the 
energy of the emitted photon and M the rest mass of the particle, more on this can be 
found in . 



When the universe is radiation dominated, this correction term is aM/n^k which 
is smaller than I/um and independent of G as well. When the universe is de Sitter 
dominated, the correction is GM j Ko?{« l/nu by hypothesis). 

The second correction term of eq. (^) is quadratic in At. It scales like the first 
correction term multiplied by At{a/a). Therefore as long as At is of the order of a few 
Compton wave lengths this second term is much smaller than the first one by virtue of 
eq. (H). By a similar analysis, one easily shows that the norm factors A±, eq. (|^) are 
given by (1 + lK,0(l)). 

Therefore, for macroscopic universe characterized by um >> 1, the non linear terms 
engendered by the change in energy carried by the operators tp can be correctly and safely 
discarded. Furthermore, the dominant correction comes from the desire of working with 
a single background characterized by the initial number of quanta um hence valid for 
both terms of the Green function. 

We now compute the corrections to the WKB approximation. We recall that we used 
this approximation in order to evaluate the kernels of gravity K{a2,ai;nM,n^), see eq. 



(|31|). Very important is the fact that Q is given in term of products of two kernels, see 
eq. (^). Indeed only a small part of the correction to the K^s will therefore affect the 
Green functions. 

As explained in the first order correction to the WKB approximation is obtained 
by replacing 7r(a, um, n^) = tTcI by 



n{a,nM,n^) = tTcH 1 + 



h 



2 r 



(51) 



As before, one has to specify the matter content of the universe to evaluate the modi- 
fication of the phase and the norm of Q induced by this correction. When the universe 
is matter dominated, the modification of phase of the Green function induced by the 
second term of eq. (|5lD is, see eq. (|^), 



( f^Planc k 1 

V M4 



1 



4 



C(l) « -MAt{a2,ai;nM,ny) — 



n 



(52) 



M , 



^In this respect, it is instructive to compare the perturbative treatments of 12 and ||T£] which are 
both appUed to a relativistic particle. Each treatment is closely related to the corresponding one in 
quantum cosmology. In [p^ , the development is made around zero momentum (this corresponds to 
isolate empty Einstein solutions) whereas in |19 the development is made around the change in the 
momentum induced by a given transition. This latter clearly corresponds to what was done in eq. (^9|). 
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where we have used eq. (p4\). This means that when eq. (§) is satisfied, the universe 
behaves so semiclassically that the part of the correction to the WKB approximation 
which modifies the Green function is irrelevant since it is smaller than the non-linear 
corrections of eq. (^) by an extra factor of I/um- 

Similarly, the change of the phase of the Green function associated with a different 
choice of the "normal ordering" of the operator d"^ in eq. (^) is equally irrelevant. 
Indeed another choice in defining would lead to an additional quantum potential 
term which would induce a change in eq. (|49|) through a change in the momentum Tia 



of the same order of the WKB corrections given in eq. (Wm)- 

In both cases, the relevant part of the corrections is given by the change of the 
correction obtained by comparing two neighbouring solutions characterized by um and 
um — 1. Indeed both the phase and the norm of Q are given in terms of the difference 
TT{a,nM, n^) — 7r(a, um — 1, n^). This is why small corrections are obtained and why we 
emphasize the necessity of the double development. As in eq. (|50D, it is the weight of 
the sources of gravity which guarantees the negligible character of these corrections. 

The third kind of corrections are associated with the spread of nAf , specified by the 
choice of the matrices Pm,Pj- To estimate these corrections and to compare them with 
the first kind just described, we write the mean quadratic spread as {AhmY = umc^m 
and similarly for An^. Then by repeating the algebra of eq. (^), one verifies that 
the change in the phase is controlled by the factor 1 + alj/nj^O^l). This result was 



obtained in JTOl. Similarly, the prefactor of Q is multiplied by 1 + iM Atau /nM- Thus for 
sufficiently small time intervals, the corrections induced by the spread are negligible for 

» 1 and ctm = Note however that the "1" in the factor riM + l which appears 

in the mean Green function defined in eq. (^) is meaningful only if MAtaM << 1. 

Had we worked with a universe described by a pure state instead of a diagonal 
density matrix, we would have obtained the additional interfering terms of eq. (^. This 
interfering term is more sensitive to gravitational corrections. Indeed, for t2 + ti — 2to 
large compare to um/ctm, the phase shifts due to the non-linearities introduce additional 
phases which will lead to the vanishing of that term upon summing over um and n^. 
Therefore, the backreaction of gravity, i.e. the dependence of the time intervals in 
the particle content of the universe (which is represented by these non linear terms), 
leads inevitably to decoherence for sufficiently long time intervals. All of this is valid 
for well peaked distribution of energy around their mean values. For wider spreads, the 
decoherence effects are even stronger. We shall return to this point after have considered 
the interactions in |^D|]. Let us just mention that the interactions among matter and 
radiation will inevitably engender spreads in proportional to um- Therefore density 
matrices which are initially too peaked will spread dynamically. Spreads in energy are 
intrinsic to cosmology. 

Finally we stress that up to now the insertion of the field operators at some radius 
has been performed from the outset as if one had an additional quantum system coupled 
to the field at that radius. However such an additional system carries some mass-energy 
which must be taken into account in the total matter hamiltonian. It will therefore 
affect gravity. To show how to treat the interactions of the internal constituents of the 
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universe is the purpose of [pO| ]. 
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6 Appendix A. 

Green functions with non- vanishing momenta 

In this Appendix we introduce an extension of the formalism in which the momenta 
of matter no longer vanish. This extention allows to compare the manner by which 
the notion of spatial displacement arises from homogeneous solutions to the manner 
we used to recover time lapses from stationary solutions. Furthermore, when dealing 
with homogeneous universes, this extention makes explicit the appearance of inertial 
coordinates as the "angles" conjugated to the conserved momenta, i.e. the "actions". 
Finally, the dynamical character of gravity, i.e. the backreaction, breaks the symmetry 
between the local representation and the momentum representation which exists when 
quantum field theory is considered in the absence of gravity. Indeed, gravity responds 
to the conserved momenta. This implies that the local representation is an approximate 
description which can only be defined a posteriori. 

Upon considering non-homogeneous gravitational and matter fields one has to deal 
with a set a four local constraints Ti.^{x) = in the place of the single scalar constraint, 
eq. (|14D- These local constraints are enforced by four Lagrange field-multipliers, A^(x), 
the local lapse and the local shifts. Up to now, these constraints have only been taken 
into account partially. For instance, in 0, the action of gravity and matter have been 
developed to quadratic order only. This leads to a hierarchy of constraints. Indeed, to 
that order in the fluctuations around an homogeneous background, the non-homogeneous 
part of the constraints 7i^(x) is linear in the fluctuations of matter and gravitational 
configurations. Instead, the homogeneous parts of the temporal constraint, H, i.e. the 
zero momentum component of 7i°(x), is quadratic in these fluctuations. The reason is 
very simple, the non- homogeneous parts of the Lagrange multipliers, N^{x), have no 
"background" contribution; only the zero momentum component of N^{x), i.e. N, has 
such. 

This has an important consequence for us. Upon dealing around an homogeneous 
background with matter fluctuations which leads to quadratic contributions in T-C^{x), 
one has to work without gravitational fluctuations if one develops the gravitational fluc- 
tuations in the action to quadratic order only. Accordingly, one should ignore completely 
the inhomogeneous part of the constraints. This quadratic approximation means that the 



19 



matter quanta do not interact directly among themselves. However they do contribute 
to the homogeneous background energy. 

We shall work in this framework. The action of our massive field is now 

Sm = IJ d'p dt [i^l - (M^ + (53) 

where the dimensionless quantity p is the conserved number of nodes of the field fluctu- 
ation. For closed three surfaces, p is an integer number. At this level, this action can 
be envisaged as describing a collection of homogeneous excitations characterized by the 
a dependent energy 



a) = + p2/a2 (54) 
Similarly the action of the massless conformally coupled field is 

S, = lJ d'k dv {{d,<j^,r - e<l>l) = \j d'k^ {a'^l - fcV^) (55) 

We can then proceed as in Section 3. The energy of gravity must annihilate the 
a-dependent total energy of the matter: 

Hg + H^= (^ ~^'^" + + n{M,p)n{p, a)+j d'k n(7, fc)^) = 

(56) 

Notice that the subtraction of the zero-point energy is now mandatory since we have an 
infinity of fluctuating modes. We emphasize that, being conserved quantities, p and k are 
gauge invariant quantities. Therefore, upon making variations of the gravitational kernel 
with respect to them, one shall define gauge invariant quantities as well. In particular, 
their conjugate variables are "angles", i.e. they do not appear in the Hamilton- Jacobi 
equation. This is how inertial coordinates are introduced and legitimized in the hamil- 
tonian framework. Furthermore, we treat n{M,p) as constants of motion as well, i.e. we 
work in the adiabatic approximation. This is valid if eq. (^ is satisfied. Therefore, the 
total wave function can be expressed, as in eq. (P^, as a sum of waves at fixed parti- 
cle number n(M, p), n(7, fc). The same is true for the kernel K{a2,ai;n{M,p),n{'~f,k)) 
defined as in eq. (|27|) . Then we can define the "WDW" Green function at fixed p in a 
universe characterized by the set of occupation numbers n{M,p'),n{'~f, k) by 

Q'^^^^'^^'''^'^'^\p,a2,ai) =fr Pn{M,p')Pn{^,k) {\a'2)ida2i^p{a2\^ [\ai)ida^i)p{ai\^ 
n(M v) 

^ -.Kia2,aun{M,p)-l,n{M,p' =^p),n{^,k)) 

I 

-.K{a2, ai; n{M, p) + 1, n(M, pV P) , k)) 



2^Vt{p,ai)Vt{p, 02 
n{M,p) + 1 



2y^Jl(p, ai)fi(p, 02) 
<— > <— > 

ida^ida^K{a2, ai; n{M,p'),n{'y, k)) (57) 
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where the derivatives act on the kernels only. 

Upon using the WKB expressions for the kernels K and to first order in the change 
in energy, one obtains 

/ 71 ^- \ « r*^ dt'n(p,a(t')) / / , ^ \ , 1 N -i dt'n(pMt')) 

yxl) [P,0-2,0-l) — , 1 / 

= G;^^'''\p,At{a,,a,;n{M,p'),n{j,k))] (58) 

where G^^^^'^\p,At{a2,ai;n{M,p'),n{'y,k))] is the Green function evaluated in a uni- 
verse whose particle content is determined by the n{M , p') , n{'j , k). 

We are now in position to determine to what extend one recovers the usual local 
representation. One first introduces x, the conjugated to the wave number p, by defining 
the (local) field operator tplx) = J Sp e'^^^ipp. Then the Green function of these operators 
is 



Pn{M,p')Pni'y,k)\a2)ida2'ip{x2) (02! |ai)z ?/'(Xi) (fli | 



It is a function of X2 — Xi only since p is a conserved quantum number. However it is not 
a Lorenz invariant function since the cosmological radius is a time dependent function. 

The validity of the local representation depends on the space-time separation of 
the two field insertions at ai,Xi and a2,X2- Indeed, the momenta which furnish the 
dominant (saddle point) contribution to g^^^'^^'^^"''''\x2,Xi,a2,ai) are centered around 
the solution of 

Ax = dp da[7r{a,n{M,p'),n{'y,k)) — 'K{a,n{M,p) — l,n{M,p' ^ p),n{'y,k))] 

J ai 

= r dt dpfiip, a) [1 + 0{p/nMM)] (60) 

Therefore the correction term is negligible only if the two points are not too close to the 
light cone defined by Ax = Ari{a2, ai; n{M,p'),n{'j, k)) and reached in the limit p —>■ oo. 
To compute the value of a massive Green function in quantum gravity on the light cone 
is meaningless since the "recoil" of gravity destroys the background. 
The lessons of this Appendix are the following 

1- There is a regime in which it is legitimate to keep only the constraint eq. (|56| ) 
and the dispersion relation eq. (0) of individual quanta. This regime correspond to 
homogeneous background driven by incoherent matter quanta, i.e. the state of the matter 
is such that the expectation value of the momentum operator Ti-lnatteri^) vanishes. 

2- In that regime, gravity is only sensitive to the various occupation numbers n{M,p) 
weighted by Q{p,a). Therefore, all matter quantities are evaluated in the momentum 
representation, at fixed p. Only then, one may introduce, by Fourrier transform, local 
representations of these quantities since the reaction of gravity is a nonlinear function 
of n{M,p) and fl{p,a). 
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3- The local representation is inevitably an approximation. Indeed, in order to local- 
ize with arbitrary precision one needs to include arbitrary momenta and thus arbitrary 
energies. But the equality between the WDW Green function ^^''^'^'^ ^'"''^''^^(p, 02, cti) 
and the corresponding Green function evaluated in the background requires that the 
exchanged energy fl{p, a) be much smaller that the background energyQ Thus the local- 
izability of an event is limited in time and in space by this restriction on the exchanged 
energy. Therefore, exactly like the fact that the time lapse is an approximate concept, it 
appears that spatial displacements X2 — Xi have a secondary approximative status. The 
distance is a useful concept which can be introduced for convenience, a maquillage to 
use the word of Banks |^. This latter aspect will become more clear upon considering, 
in the next paper the interactions of an heavy particle and the recoils induced by 



emissions of light quanta. Indeed, without interactions, there is no physical justification 
to introduce the local operator ipi^x) = J d^p e^^^ipp. 

7 Appendix B. The corrections in v/V in a micro- 
canonical ensemble 

In the Appendix we shall compute the induced (almost canonical) partition function of a 
small system pertained in a bigger one. The total system is microcanonically distributed. 
For simplicity, we shall suppose that the both systems are homogeneous and that the 
contact energy can be neglected. The total system has a volume V + v and the little 
system has a volume v. 

The goal of this exercise is to determine how to obtain a well defined expansion in 
v/V. This is not trivial and bears many resemblance with the debated extraction of the 
background in quantum cosmology. 

The total system has an energy E. Thus the number of states (to be more precise, 
the density of states) can be expressed as 

^totaiiE) = r de Qv{E - e) ^(e) (61) 
Jo 

where Qv{E — e) is the number of states when the rest of the big system with volume 
V has an energy E — e and where u^e) is the number of states of the little system with 
energy e. 

Two alternative ways to evaluate this quantity will be considered and compared. 
The first one consists in neglecting completely the energy of the little system (this seems 
legitimate since the ratio of the mean energies e/ (E—e) scales like v/V) and in developing 
around the vacuum solution e = 0. To first order in the change of flv{E — e) in e, one 
obtains 

QtotamiE) - Qv{E) de e-^^'u{e) 

Jo 

= nv{E)z{P,) (62) 



^ Moreover upon considering the set of local constraint, this restriction on H, might be re-enforced. 
We hope to return on that point in a next publication 
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where the inverse temperature (3i is defined by (3i = —d^ln{Qv{E — e)) evaluated at 
e = and where z{l3) designates the canonical partition function of the small system. 

The second method consists in taking into account the the mean energy of the little 
system e = e determined by the saddle point of eq. (|6lD. This saddle expresses, as usual, 
the equality of the temperatures: 

dJnnviE - e) + d,lnuj{e) = (63) 

where 

/32 = dE\nnviE-e) (64) 

The derivative is now evaluated at the saddle point energy e = e. Then, to first order in 
the change in energy, i.e. e — e, the partition function for the little system reads 



Jo 

= nv{E-e) e^^' z{p2) (65) 

To this linear order in the development of the density of the big system f^y, the 
only difference between the two approaches lies in the difference in the temperatures 
only. Furthermore, one finds that /?i — ^2 scales like v/V. Therefore if one takes the 
thermodynamical limit at this moment, i.e. E — > oo, — > oo with E/V fixed, both 
methods give the same results. 

However, this is no longer the case if one wants to evaluate finite size effects, that is 
the fact that the big system has a finite heat capacity. Then, the difference between the 
two developments show up to second order in the energy fluctuation. (This is inevitable 
since /3i — P2 is of the order of v/V.) Indeed, by developing Qv{E — e) to second order 
in e, one has 

Qtotai-AE) = ^v{E) r de e'P^' e-^^'/2^-a;(e) (66) 
Jo 

in the flrst case and 

^totai;2{E) =nv{E- e) de e-^^(^-^") e-^'^'-'^'/'^'-uie) (67) 

Jo 

in the second case. We have neglected the difference in the heat capacity of the big system 
(deflned by d^\nfl{E — e) = —(3'^/Cy) which is induced by the different "background" 
energies E and E — e, since Cy appears in the correction term only. 

It is convenient to introduce the entropies s = lna;(e), S = \iaQ{E — e) as well as 
their relations with the mean energies: s = ^i[3e = C,2vP^ and S = ii[3{E — e) = ^2^^/?^- 
In those equations, the factors are dimensionless and of the order of 1. It is also 
convenient to relate the speciflc heats to the entropies of both systems as Cy = ^38 and 
Cy = = ((e — e)^)/5^. Then one flnds that the mean correction in the flrst case scales 
like 
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while in the second case it scales only like 



These results follows from the fact that, while the temperature is essentially fixed by 
the big system, it is the little one which determines, at a given temperature, what are 
the mean quadratic fluctuations, i.e. the width around the saddle energy e. Notice that 
there is no need to differentiate the temperatures to that order in v/V . 



By inspection of eqs. (^, |6^), one sees that only the second development is really 
a power series in v/V. Instead, in the first case, the corrections contains the factor 
v'^(3'^ /V which is ambiguous because it can be much bigger or much smaller than one 
when V « V. Indeed, it suffices to consider a system such that s « S but with 
> S. 

The lesson is that when one wants to compute finite size effects, i.e. first order 
correction in v/V, it is mandatory to take the mean energy of the little system into 
account upon evaluating the temperature, the zeroth order "classical" solution which 
corresponds to the background from which the temperature is computed. The paral- 
lelism between these developments and the one appearing in eqs. (^, ^9|) is manifest. 
*corresponds zeroth order "classical" solution which In this respect we want to point 
out that the factor e^^^ of eq. (^) finds its perfect counterpart in ref. ||l^ wherein the 



Schroedinger time dependent equation differs from the usual one by e**"^. Its origin is 
identical to our, i.e. the background (equilibrium) is determined through a saddle point 
which does involve the mean "matter" energy e. 
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